Submitted to the Astrophysical Journal, January 10, 2012 

Preprint typeset using L^T^X style cmulatcapj v. 5/2/11 



CLASH: PRECISE NEW CONSTRAINTS ON THE MASS PROFILE OF ABELL 2261 

Dan Coe 1 , Keiichi Umetsu 2 , Adi Zitrin 3 , Megan Donahue 4 , Elinor Medezinski 5 , Marc Postman 1 , Mauricio 
Carrasco 6,7 , Timo Anguita 6,8 , Margaret J. Geller 9 , Kenneth J. Rines 910 , Antonaldo Diaferio 1112 , Michael 
J. Kurtz 9 , Larry Bradley 1 , Anton Koekemoer 1 , Wei Zheng 5 , Mario Nonino 13 , Alberto Molino 14 , Andisheh 
Mahdavi 15 , Doron Lemze 5 , Leopoldo Infante 6 , Sara Ogaz 1 , Peter Melchior"', Ole Host 17 , Holland Ford 5 , 
Claudio Grillo 18 , Piero Rosati 7 , Yolanda Jimenez-Teja 14 , John Moustakas 19 , Tom Broadhurst 20 Begona 
Ascaso 14 , Ofer Lahav 17 , Matthias Bartelmann 3 , Narciso Beni'tez 14 , Rychard Bouwens 21 , Or Graur 2 , Genevieve 
Graves 23 , Saurabh Jha 24 , Stephanie Jouvel 17 , Daniel Kelson 25 , Leonidas Moustakas 26 , Dan Maoz 22 , Massimo 
Meneghetti 2 ' , Julian Merten 25 , Adam Riess 1,5 , Steve Rodney 5 , and Stella Seitz 28 

Submitted to the Astrophysical Journal, January 10, 2012 

ABSTRACT 

We precisely constrain the inner mass profile of Abell 2261 (z = 0.225) for the first time and 
determine this cluster is not "over-concentrated" as found previously, implying a formation time in 
agreement with ACDM expectations. These results are based on strong lensing analyses of new 
16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble 
(CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5- 
band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M V i r = 
2.2 ±0.2 x l0 15 M(7)hyQ (within r OT > ~ 3 Mpc h^) and concentration c V i r = 6.2 ±0.3 when assuming a 
spherical halo. This agrees broadly with average c(M, z) predictions from recent ACDM simulations 
which span 5 < (c) < 8. Our most significant systematic uncertainty is halo elongation along the line 
of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations 
and find it to be ~ 35% lower than our lensing-derived profile at r25oo ~ 600 kpc. Agreement can be 
achieved by a halo elongated with a ~ 2 : 1 axis ratio along our line of sight. For this elongated halo 
model, we find M V i r = 1.7 ± 0.2 x 10 15 Mgh^ and cw = 4.6 ± 0.2, placing rough lower limits on 
these values. The need for halo elongation can be partially obviated by non-thermal pressure support 
and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of 
background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower 
M vir further by ~ 7% and increase c V i r by ~ 5%. 

Subject headings: cosmology: dark matter — galaxies: clusters: individual (Abell 2261) — gravita- 
tional lensing: strong — gravitational lensing: weak — cosmology: dark energy — 
galaxies: evolution 
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1. INTRODUCTION 

Detailed observational constraints of dark matter halos 
yield important tests to our understanding of structure 
formation (Natarajan et al. 2007; Borgani & Kravtsov 
2011), the particle nature of dark matter (Clowe et al. 
2006; Kceton & Moustakas 2009), and perhaps the na- 
ture of dark energy as well (Grossi & Springel 2009). 
Large cluster surveys require precisely determined clus- 
ter masses to calibrate their observables and achieve their 
full potential to constrain cosmology (Henry et al. 2009; 
Allen et al. 2011). 

The galaxy clusters studied best via gravitational lens- 
ing appear to have more densely concentrated cores 
than clusters of similar mass and redshift formed in 
ACDM simulations (Broadhurst et al. 2008; Broadhurst 
k Barkana 2008; Oguri et al. 2009; Richard et al. 2010; 
Sereno et al. 2010; Zitrin et al. 2011a; Postman et al. 
2011). (See results from other methods reviewed in Fcdcli 
2011 and Bhattacharya et al. 2011.) Some of this dis- 
crepancy is due to bias, as the clusters selected for the 
most detailed lensing studies were among the strongest 
gravitational lenses known. However it is estimated that 
even this large (~ 50 - 100%) bias cannot fully explain 
the high observed concentrations (Hcnnawi et al. 2007; 
Oguri et al. 2009; Meneghetti et al. 2010, 2011, although 
see Oguri et al. 2011). Baryons, absent from these dark 
matter only simulations, are found to only modify clus- 
ter concentrations at the < 10% level (Duffy et al. 2010; 
Mead et al. 2010; Fedcli 2011). 

If confirmed, this result would imply that galaxy clus- 
ters formed earlier than their counterparts in simulated 
ACDM universes. We expect that the higher density of 
the earlier universe would remain imprinted on the clus- 
ter cores as we observe them today (e.g., Jing 2000; Bul- 
lock et al. 2001; Wechsler et al. 2002; Zhao et al. 2003). 

Another possible hint of early cluster formation may 
be galaxy clusters detected at z > 1 which are perhaps 
unexpectedly massive (Stanford et al. 2006; Eisenhardt 
et al. 2008; Jee et al. 2009; Huang et al. 2009; Rosati et al. 
2009; Papovich et al. 2010; Schwope et al. 2010; Gobat 
et al. 2011; Jee et al. 2011; Foley et al. 2011; Santos 
et al. 2011a, b; Planck Collaboration et al. 2011). How- 
ever one should be cautious about the statistical interpre- 
tation of such results (Paranjapc et al. 2011; Hotchkiss 
2011; Hoyle et al. 2011; Waizmann et al. 2011; Harri- 
son <fc Coles 2011a, b). Building on results from large X- 
ray surveys (Vikhlinin et al. 2009; Mantz et al. 2010b), 
large new Sunyaev-Zel'dovich (SZ) surveys will continue 
to constrain cosmology based on cluster abundance mea- 
surements as functions of mass and redshift (Sehgal et al. 
2011; Benson et al. 2011). 

Mechanisms proposed to explain such early growth in- 
clude departures from an initially Gaussian spectrum of 
density fluctuations (e.g., Chongchitnan & Silk 2011), 
though we note some such non-Gaussian models can be 
ruled out based on cosmic X-ray background measure- 
ments (Lemze et al. 2009). Early growth may also be 
explained by higher levels of dark energy in the past. 
This idea, dubbed Early Dark Energy (EDE; Fedeli & 
Bartelmann 2007; Sadeh & Rephaeli 2008; Francis et al. 
2009; Grossi & Springel 2009), would have suppressed 
structure growth in the early universe, such that clus- 
ters would have had to start forming sooner to yield the 



numbers we observe today. Other dark energy theories 
with similar implications have also been proposed (e.g., 
Baldi 2011; Carlesi et al. 2011). 

Significant improvements in these observational con- 
straints are being obtained by CLASH, the Cluster Lens- 
ing and Supernova survey with Hubble (Postman et al. 
2011). 29 CLASH is a 524-orbit multi-cycle treasury HST 
program to observe 25 galaxy clusters (0.18 < z < 0.89) 
each in 16 filters with the Wide Field Camera 3 ( WFC3; 
Kimble et al. 2008) and the Advanced Camera for Sur- 
veys (ACS; Ford et al. 2003) over the course of three years 
(HST cycles 18-20). Importantly, 20 CLASH clusters 
were X-ray selected to be massive and relatively relaxed. 
This avoids the strong bias toward high concentrations in 
previously well-studied clusters selected for their lensing 
strength. 

Abell 2261 (hereafter, A2261) was observed as part of 
the CLASH program. It has a redshift of z = 0.2249 as 
measured by Crawford et al. (1995) and refined by Rincs 
et al. (2010). 

Weak lensing (WL) analyses of ground-based imaging 
of A2261 (Umctsu et al. 2009; Okabe et al. 2010) yielded 
concentration measurements of c V i r ~ 6 or ~ 10, with the 
broad range attributed to measurement uncertainties, 
the details of the analysis method used, and perhaps sub- 
ject to uncertainty due to massive background structures 
identified at z ~ 0.5. The latter value (c V i r ~ 10) would 
be significantly higher than predicted for an average re- 
laxed cluster of A2261's mass and redshift: c v i r ~ 5 from 
Duffy et al. (2008), although an analysis of more recent 
simulations yields a much higher prediction: c V i r ^8.5 
(Prada et al. 2011). The WL measurements had over- 
lapping uncertainties, but a preliminary strong lensing 
(SL) measurement of the Einstein radius (Re ~ 40" for 
a background source at z s = 1.5) supported the larger 
value with smaller uncertainties: c V i r — 11 ± 2 (Umetsu 
et al. 2009). This result was also included in Oguri et al. 
(2009) as one of ten well-studied clusters, all of which 
had higher than predicted concentrations. 

In this work, we revisit both the strong and weak 
lens modeling. Our deep 16-band HST imaging reveals 
strongly lensed (multiply imaged) galaxies all undetected 
in the previous HST imaging (0.5-orbit WFPC2 F606W) 
and allows us to derive robust and precise photometric 
redshifts for these arcs, a key ingredient for our mass 
model of the cluster core. 

Detailed strong lensing analysis is required to accu- 
rately and precisely measure the inner mass profile and 
concentration of A2261. By probing the mass profile 
over a combined two decades of radius, joint analysis of 
strong plus weak lensing yields significantly higher preci- 
sion measurements of cluster virial masses and concentra- 
tions than either method alone (Meneghetti et al. 2010). 

This paper is organized as follows. We describe our 
HST (§2) and MMT spectroscopic (§3) observations fol- 
lowed by our strong lens mass modeling (§4). We then 
introduce our ground-based imaging and weak lensing 
analyses (§5) and derive joint strong + weak lensing con- 
straints (§6). We constrain halo triaxiality from joint 
lensing + X-ray constraints in §7 and finally compare 
our mass profile with results from simulations in §8. The 
formation time of A2261 is discussed in a broader con- 

29 http:/ /www. stsci.edu/~postman/CLASH 
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Table 1 

CLASH HST Observations of the Core of A2261 



Camera / 


Filter 


HST 


Exposure Time 


Channel 


Element 


Orbits 


(sec) 


Wr to/ U Vlb 


r 220 W 


1.5 


3671 


TT7T7i/^io /TT"\7'TG 

W-b Co/ U Vlb 


b 275 W 


1.5 


3745 


\XTT?f^') /TT'\7'TG 

WrLo/ U Vlb 


r 3odW 


1.0 


2408 


W r \_jq/ U V lb 


r oyUW 


l.U 


Z40D 


ACS/WFC 


F435W 


1.0 


2077 


ACS/WFC 


F475W 


1.0 


2064 


ACS/WFC 


F606W 


1.0 


2057 


ACS/WFC 


F625W 


1.0 


2064 


ACS/WFC 


F775W 


1.0 


2072 


ACS/WFC 


F814W 


2.0 


4099 


ACS/WFC 


F850LP 


2.0 


4148 


WFC3/IR 


F105W 


1.0 


2814 


WFC3/IR 


F110W 


1.0 


2514 


WFC3/IR 


F125W 


1.0 


2514 


WFC3/IR 


F140W 


1.0 


2411 


WFC3/IR 


F160W 


2.0 


5029 


Note. 


Parallel observations 


and supernova 



follow-up observations are not described here nor uti- 
lized in this work. 

text including other observational probes in §9, and we 
summarize our conclusions in §10. 

Where necessary to calculate distances, etc., we assume 
a concordance ACDM cosmology with h = 0.7, f2 m = 
0.3, n A = 0.7, where H = 100 h km s^ 1 Mpc" 1 . In this 
cosmology, at A2261's redshift of z = 0.225, 1" w 3.59 
kpc h^Q Rj 2.51 kpc where h = O.7/170. Furthermore 
at this redshift, the cluster virial radius is defined as that 
which contains an average overdensity of A c « 115 times 
critical, where A c w 187T 2 — 82Ha(z) — 39fl\ (z) based on 
spherical collapse theory (Bryan & Norman 1998). 

2. HST OBSERVATIONS 

We observed A2261 (BCG R.A. = 17 h 22 m 27:2, deck 
= +32°07'57" [J2000]) as part of the CLASH program in 
HST Cycle 18 between 2011 Mar 9 and May 21 to a total 
depth of 20 orbits in 16 WFC3 and ACS filters, spanning 
~ 2,000A - 17,000A (Table 5; GO 12066; PI Postman). 
The images were processed for debias, flats, superflats, 
and darks using standard techniques, and then co-aligned 
and combined using drizzle algorithms. See Koekemoer 
et al. (2007) and Postman et al. (2011) for details. 

In order to better reveal faint lensed images, we mod- 
eled and subtracted the BCG light in all 12 ACS+IR fil- 
ters. We used the isophote fitting routine, SNUC, which 
is part of the XVISTA image processing system, to derive 
two-dimensional models of the bright early-type galaxies 
in A2261, including the BCG. SNUC is capable of simul- 
taneously obtaining the best non-linear least-squares fits 
to the two-dimensional surface brightness distributions in 
multiple, overlapping galaxies (Laucr 1986). The models 
were derived independently for each CLASH passband. 
Fits were performed using concentric isophotes but the 
position angles and the ellipticities of the isophotes were 
left as free parameters. The models were then subtracted 
from the original image to produce a bright-galaxy sub- 
tracted image. 

We used SExtractor (Bertin & Arnouts 1996) to detect 
objects and measure their photometry. For arcs elud- 
ing this initial detection, we constructed manual aper- 
tures which were then forced back into SExtractor using 
SExSeg (Coe et al. 2006). Isophotal apertures were used 



as they have been shown to yield robust colors (Bem'tez 
et al. 2004). 

Based on this photometry, we measured photometric 
redshifts using BPZ (Bem'tez 2000; Bem'tez et al. 2004; 
Coe et al. 2006). Spectral energy distribution (SED) 
templates are redshifted and fit to the observed photom- 
etry. A Bayesian analysis tempers the qualities of fit with 
a prior: the empirical likelihood of redshift as a function 
of both galaxy magnitude and type (e.g., bright and/or 
elliptical galaxies are rare at high redshift). Here we 
used 11 SED templates originally from PEGASE (Fioc & 
Rocca-Volmerange 1997) but strongly recalibrated based 
on photometry and spectroscopic redshifts of galaxies in 
the FIREWORKS catalog Wuyts et al. (2008). These 
templates yield < 1% outliers for high quality spectro- 
scopic samples and therefore implicitly encompass the 
full range of metallicities, extinctions, and star forma- 
tion histories of real galaxies. 

The color images used in this paper were produced 
automatically using the publicly available Trilogy soft- 
ware. 30 Trilogy determines the intensity scaling auto- 
matically and independently in each color channel to 
display faint features without saturating bright features. 
The scalings are determined based on a sample of the 
summed images and two input parameters. One sets the 
output luminosity of "the noise", currently determined 
as 1-er above the sigma-clipped mean. The other pa- 
rameter sets what fraction of the data (if any) in the 
sample region should be allowed to saturate. Default 
values for these parameters (0.15 and 0.001%, respec- 
tively) work well, but the user is able to tweak them. 
The scaling is accomplished using the logarithmic func- 
tion y — a\og w (kx + 1) clipped between and 1, where 
a and k are constants determined based on the data and 
desired scaling parameters as described above. 

3. MMT/HECTOSPEC SPECTROSCOPY 

The Hectospec instrument mounted on the 6.5-meter 
MMT is a multiobject fiber-fed spectrograph with 300 
fibers deployable over a 1° diameter field (Fabricant et al. 
2005). We targeted probable A2261 cluster members 
based on their proximity to the expected cluster red se- 
quence and proximity to the BCG. In addition, all targets 
were sufficiently bright with magnitudes 16 < r' < 21 in 
SDSS DR7 (Abazajian et al. 2009). To observe all targets 
we used four pointings each with 3x20 minute exposures 
within the Hectospec queue schedule. 

After processing and reducing the spectra, we used 
the Hectospec pipeline (Mink et al. 2007) based on the 
IRAF package RVSAO (Kurtz & Mink 1998) to cross- 
correlate the spectra with a set of standard Hectospec 
templates. In total, we obtained high quality redshifts 
for 572 galaxies, including 308 within our Subaru analy- 
sis region (within 17.8' of the BCG). These redshifts are 
provided in Table 3. Scattered light from bright stars in 
this field limits spectroscopy of this system. 

Redshift uncertainties are measured for each galaxy in- 
dividually and have been empirically quantified globally 
as follows. The SHELS survey (Geller et al. 2005) carried 
out with Hectospec obtained repeat observations of 1468 
galaxies, yielding a mean internal error of 56 km s _1 for 
absorption-line objects and 21 km s _1 for emission-line 

30 http:/ /www. stsci.edu/~dcoc/trilogy/ 
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Figure 1. Multiple images of background galaxies strongly lensed by A2261. All were identified in this work based on our deep, 
multiband HST imaging and lens modeling. Each is located directly above and to the left of its label. Close ups may be seen in Fig. 2. 
The prominent arc marked "X" is not multiply-imaged according to our models. Lensing critical curves from our primary lens model are 
plotted for background sources at redshifts z B = 1 (blue), 2 (green), and 7 (red). These HST ACS+WFC3/IR color images were produced 
automatically using Trilogy (§2) which reveals faint features without saturating bright areas such as the BCG core. Filters were assigned 
colors as follows: Blue = F435W + F475W; Green = remaining ACS filters; Red = WFC3/IR. The green "figure 8" patterns are ACS 
reflection artifacts from a bright star to the SE. North is up, East is left. Left: The diamond-shaped hole in the ACS images was filled in 
with an IR image tinted yellow to roughly match the color of this portion of the BCG's stellar halo. Right: Light has been modeled and 
subtracted from the BCG and a few other cluster galaxies close to the arcs (§2). The residual pattern near the location of the BCG center 
(aside from the bright knots) is most likely due to a combination of small model artifacts and real asymmetries in the stellar distribution. 




Figure 2. Multiple images of galaxies strongly lensed by A2261, as identified in this work (see §4.2, Figure 1, and Table 2). The BCG 
has been modeled and subtracted from all of these images. Most are ACS color images (B = F435W+F475W, G = F606W+F625W, R = 
F775W+F814W+F850LP), while those at bottom right use WFC3/IR filters (B = F105W+F110W, G = F125W+F140W, R = F160W). 
Each image stamp is 5.2" 19 kpc on a side. 



CLASH: Precise Mass Profile of Abcll 2261 



■5 



Table 2 

Multiple Images of Galaxies Strongly Lensed by Abell 2261 



ID 


R.A. 
(J2000.0) 


Decl. 
(J2000.0) 


Magnitude 
F775W (AB mag) 


Photometric 
Redshift" 


Lens Model 
Redshift 6 


la 
b 


17 22 25.452 
17 22 27.295 


+32 08 25.02 
+32 07 42.58 


23.02 + 0.04 
22.82 + 0.04 


* 4 qq+0-03 
n 4s+ - 03 

°- 48 -0.03 


4.44 


2a 
b 

c 


17 22 26.338 
17 22 26.133 
17 22 25.897 


+32 08 20.54 
+32 08 18.59 
+32 08 17.16 


26.81 ± 0.14 
25.78 + 0.16 
27.04 + 0.20 


*3.89t^i3 
4 07+° 11 
o 74+0.14 
J - '^-0.13 


3.81 

11 


3a 
b 
c 
d 


17 22 28.917 
17 22 28.912 
17 22 26.972 
17 22 25.391 


+32 07 55.13 
+32 07 53.96 
+32 07 38.49 
+32 08 18.13 


24.69 + 0.07 

11 C 

26.63 + 0.28 

27.64 + 0.28 


26 +2,ss 

55 C 

2.60 + °-^ 
*3 S8+ ' 11 


3.24 

11 
11 


4a 
b 

c 


17 22 28.569 
17 22 28.369 
17 22 25.217 


+32 08 08.00 
+32 07 37.51 
+32 08 13.52 


24.63 ± 0.04 
24.12 + 0.04 
24.17 + 0.05 


„ .„-!-() OS 

3-48t ; ^ 

o no+0.07 

" J -" JO -0.07 
*^ 40+ 


3.31 

11 


5a 
b 

c 


17 22 28.589 
17 22 

17 22 25.202 


+32 08 07.15 

1 OO /""IT OT OO 

+32 07 37.38 
+32 08 13.71 


25.23 + 0.09 

o tr r\ a i n no 

25.04 ± 0.08 
24.91 ± 0.09 


°- a -0.25 
°- 29 -0.04 


3.31 

11 


6a 
b 
c 


17 22 27.069 
17 22 27.115 
17 22 28.615 


+32 08 09.30 
+32 08 01.76 
+32 07 28.74 


24.49 + 0.07 
21.79 + 0.01 
24.58 + 0.07 


3.26±r 07 
22+ - 01 

*o. 94+0.03 
-0.03 


3.09 

11 


7a 
b 


17 22 26.599 
17 22 29.326 


+32 07 51.81 
+32 07 59.35 


24.31 ± 0.07 
23.06 + 0.05 


i c 7 +0.01 
1 -°'-0.02 
*i =4+0. 01 


1.74 

11 


8a 
b 


17 22 26.158 
17 22 26.737 


+32 08 25.87 
+32 07 41.61 


25.87 + 0.14 
27.17 + 0.17 


*4 q9+ ula 
-0.13 

4 82+ ' 05 


4.93 

11 


9a 
b 


17 22 25.304 
17 22 27.223 


+32 08 21.19 
+32 07 39.20 


27.11 ± 0.37 
26.41 ± 0.19 


*o.78iS:iS 

4.67t°; 1 i 3 5 


4.46 

11 


10a 
b 


17 22 26.379 
17 22 26.261 


+32 08 01.89 
+32 08 01.88 


24.01 ± 0.05 
24.29 + 0.07 


*1 70+0.01 

L - ' y -0.23 


1.67 


11a 
b 


17 22 26.854 
17 22 26.875 


+32 08 06.50 
+32 08 05.01 


22.43 + 0.03 

55 c 


*•? ss+ u - 04 
3.88_ 0>04 

55 c 


3.90 

11 


12a 
b 


17 22 26.563 
17 22 29.045 


+32 08 12.81 
+32 07 34.33 


26.42 + 0.21 
27.50 ±0.15 b 


9 A 9 + 0. Ob 
z -^ z -0.70 
* 9 4 o+0.34 
z -^°-0.26 


2.83 

11 



a The best isolated, least contaminated arcs (marked with *) are chosen to provide the 
input rcdshift for each system. These arc then optimized by the lens model. Note robust 
photometric redshifts are not expected for arcs fainter than ~26th magnitude. Uncertain- 
ties are formally 95% C.L. Broad, asymmetric error bars generally indicate a bi-modal 
solution. 

b Results from LensPerfect 

c Continuous arcs were analyzed with a single photometric aperture for more robust pho- 
tometric rcdshift estimates. 
d F814W 
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Figure 3. Spectral energy distributions (SEDs) of arcs as observed in our 16 HST filters. Also plotted in gray are BPZ SED fits to the 
best isolated arcs marked with asterisks (*) both in this figure and in Table 2. For each arc, the SEDs of all images should have similar 
shapes though these may shift vertically in magnitude as their magnifications may vary. The photometry of a few faint images may be 
contaminated by cluster light, despite our best efforts to model and subtract the BCG and other cluster galaxies. Most notably, 6b is a 
faint image very near the BCG core. 




Figure 4. Two strong lensing mass models derived for A2261 exhibiting our range of density profile slopes (plotted in Fig. 5). The steepest 
of five models from the Zitrin et al. (2009b) method is shown at left and a shallower model from LensPerfect (Coe et al. 2008) is shown at 
right. The former assumes light approximately traces mass while the latter does not except for a prior on the peak density position. The 
LensPerfect model is only well constrained within the white polygon which traces the outermost multiple images. Plotted are contours of 
projected mass density in units of the strong lensing critical density (ft = 1) for a background object at z s = 2 (~ 1.9 X 1O 15 M0/Mpc 2 ). 
The contours are logarithmically spaced such that each contour is 10% denser than the next. Both ACS+IR color images are to the same 
scale: 78" ~ 280 kpc on a side. 
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Table 3 

MMT/Hectospec Spectroscopic Redshifts for Galaxies within the 
Subaru FOV 



R.A. 


Decl. 




Redshift 


Cluster 


(J2000 deg) 


(J2000 deg) 


Redshift 


Uncertainty 


Member" 


260.29303 


32.07600 


0.22096 


0.00017 





260.30228 


32.26368 


0.38051 


0.00020 





260.32202 


32.27190 


0.13698 


0.00007 





260.32373 


32.22027 


0.11367 


0.00009 





260.34888 


32.04460 


0.15224 


0.00008 






Note. — This table is published in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 
a Based on dynamical analysis of the cluster caustics. 



objects (Fabricant et al. 2005). Comparison of spectro- 
scopic redshifts obtained for 379 galaxies in both SHELS 
and SDSS DR7 yields Av = 10 km s" 1 ± 35 km s" 1 . 
Note that Av = 30 km s _1 corresponds to Az sa 0.0001. 

Based on the measured redshifts, we identified clus- 
ter members using the caustic technique (e.g., Diaferio 

2009) . The technique locates two curves, the caustics, in 
the cluster redshift diagram, namely the plane of the linc- 
of-sight velocities of the galaxies versus their projected 
clustercentric distances. The caustics are related to the 
escape velocity from the cluster and define an area of 
the redshift diagram where most of the cluster members 
reside. Samples of members identified with the caustic 
technique are at least 95 percent complete and contam- 
inated by interlopers by 10 percent at most (Serra et 
al. 2012, in preparation). This procedure has been used 
to identify galaxy members of clusters and groups (e.g., 
Rines et al. 2005; Rines & Diaferio 2006, 2010) as well 
as stellar members of the Milky Way halo (Brown et al. 

2010) and of dwarf spheroidals (Serra et al. 2010). 

We applied the technique to a sample of 777 spectro- 
scopic redshifts in the field of A2261 (Rines et al. 2012, 
in preparation). The technique identifies 209 members 
within 6.6 Mpc h^ of the cluster center. We note that 
the caustic technique also locates the cluster center ~ 6' 
south of the BCG location and at a redshift 0.0017 larger. 
This result indicates that the dynamical structure of 
A2261 might be more complex than expected. 

We used cluster members identified here to validate 
galaxy selections based on broad-band photometric col- 
ors used in our strong (§4.1) and weak lensing analyses 
(§5.5). 

4. STRONG LENS MASS MODELING OF THE CLUSTER 

CORE 

We performed three semi-independent strong lensing 
(SL) analyses on the A2261 HST images. We used the 
method of Zitrin et al. (2009b) to perform the primary 
strong lensing analysis (§4.1) including the identification 
of multiple-image systems (§4.2). These multiple images 
are all identified for the first time in this work based on 
our HST imaging and lens modeling. We verified these 
identifications using the Lenstool 31 modeling software 
(Kneib 1993; Julio et al. 2007) as well as LensPerfect 32 
(Coe et al. 2008, 2010), a "non-parametric" method that 
does not require the assumption that light traces mass 

31 http://www.oamp.fr/cosmology/lenstool/ 

32 http:/ /www. its. caltech.edu/~coe/LcnsPerfect/ 



(§4.3). Finally we combined the results from all methods 
yielding an average cluster core mass profile with uncer- 
tainties (§4.4). By utilizing various modeling methods, 
we captured the true systematic uncertainties more reli- 
ably than generally possible using a single method. 

4.1. Primary Strong Lensing Analysis Method 

The Zitrin et al. (2009b) method was adapted from 
that used in Broadhurst et al. (2005a), reducing the 
number of free parameters to six, and has been used 
extensively since (Zitrin et al. 2009a,b; Zitrin & Broad- 
hurst 2009; Zitrin ct al. 2010, 2011a,b,c,d,e; Merten et al. 
2011). Basically, the mass model consists of three com- 
ponents: the cluster galaxies, a dark matter halo, and 
an external shear to account for additional ellipticity in 
the mass distribution in the plane of the sky within or 
around the core. Cluster galaxy light is assumed to ap- 
proximately trace the dark matter; the latter is modeled 
as a smoothed version of the former, as described below. 

We identified 118 probable cluster galaxies along the 
"red sequence" which is well isolated in F814W-F475W 
color-magnitude space. We verified, using additional fil- 
ters and photometric redshifts, that this selection is ro- 
bust. We also compared this selection to Hectospec spec- 
troscopic redshifts available for 15 galaxies within the 
HST FOV. We correctly identified 11 of the 13 cluster 
members, missing one near the FOV edge and another 
near the bright star. We incorrectly identified one of the 
two foreground objects (z — 0.1693) as a cluster member 
as it fell along our red sequence. These three particu- 
lar misidentifications have a negligible effect on our mass 
model as they all lie at R > 80", well outside the strong 
lensing region where multiple images are formed. Nor do 
these ~ 10% rates of incompleteness and contamination 
significantly affect our mass profile as evidenced in part 
by our other analyses (§4.4). The cluster members pro- 
vide a parameterization for the mass model which is not 
required to be exact but rather provides a starting point 
which is molded to fit the data. 

Each cluster galaxy is modeled as a power law density 
profile, its mass scaling with flux observed in F814W. 
This mass distribution is then smoothed using a 2D poly- 
nomial spline to provide a model for the dark matter 
distribution in the cluster halo. This "smooth" mass 
component is added to the more "lumpy" (unsmoothed) 
galaxy component. Finally, an external shear is added. 

In all, there are 6 free parameters: the mass sear- 
ings of both the "smooth" and galaxy components, the 
power law of the galaxy density profiles, the degree of the 
smoothing polynomial, and the amplitude and direction 
of the external shear. The routine iterates over lens mod- 
els to find that which best reproduces the observed po- 
sitions of the strongly lensed images in the image plane, 
rather than in the source plane, which can bias solutions 
toward flatter profiles and higher magnifications. Com- 
plete details may be found in Zitrin et al. (2009b). 

In this paper, we introduce an alternative Gaussian 
convolution kernel to produce the "smooth" mass com- 
ponent. In this case, the free parameter is the Gaussian 
width instead of the polynomial spline degree. 

4.2. Multiple Images of Strongly Lensed Galaxies 

Using this method, we identified 30 multiple images 
of 12 background galaxies strongly lensed by A2261 (see 
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Figures 1 and 2 and Table 2). These are all identified for 
the first time in this work. 

We used an iterative process to identify images and add 
them to the model, beginning with those which are most 
confident. Our most confident multiple image system is 
the "claw" or U-shaped object, system 1. The distinctive 
morphology is apparent in both images, including a color 
gradient best viewed in the IR color images with the 
BCG subtracted (see Fig. 2). Image la yields a photo-z 
z ~ 4.4. The IR flux of image lb appears to be biased a 
bit high by contaminating light from the BCG (Fig. 3), 
such that the best fit SED is an early type galaxy at 
z ~ 0.5. The irregular morphology is not consistent with 
an early type galaxy. 

Based on this system, we obtained the initial mass 
model, enabling us to predict the lensed positions of 
countcrimages of other galaxies by delensing them to 
their putative true source positions and then relensing 
them with our model. Candidate counterimages were 
identified as being near the observed position, with the 
predicted lensed morphology and orientation as well as 
consistent colors and photometric redshifts. Observed 
photometry and SED fits are shown in Fig. 3. Our mul- 
tiple images generally have consistent observed SEDs (al- 
lowing for variations in magnification) and thus photo- 
z's. However some images yield unreliable photo-z's if 
they are faint and/or their light is contaminated by a 
bright nearby cluster member. To date, no spectroscopic 
redshifts are available for these galaxies strongly lensed 
by A2261. 

Our 6-parameter model is fully constrained by the posi- 
tions of our 30 multiple images. We produced four mass 
models spanning the range of profile slopes allowed by 
the data as dictated by the density power law of the 
cluster galaxies (1.1 < q < 1.2, where the surface den- 
sity £ oc R~ q ) and the polynomial smoothing degree 
(4 < S < 8) . We also generated a fifth mass model using 
a Gaussian convolution kernel with an optimized width 
of 9.1" and q = 1.3. This model is shown in Fig. 4. 

We note that our lens models do not predict counter- 
images for the large prominent arc marked with an "X" 
in Fig. 1. Instead they predict a single highly-distorted 
image, as observed. We measure its photometric redshift 
to be z = 1.19±g;g| (95% C.L.). If it were at a slightly 
higher redshift z ~ 1.5, some of our models would pre- 
dict a radial arc counterimage on the opposite side of 
the BCG core. We detect no such image. This search is 
aided by the fact that the arc is significantly detected in 
F390W where most of the other arcs drop out and the 
BCG light is significantly reduced. 

We initially identified a possible counterimage to this 
large arc with similar colors and photo-z at RA, Dec 
(J2000) = 17:22:29.4, +32:07:29 (near the bottom left 
corner of Fig. 1). However, the lens models required to 
reproduce this counterimage were significantly stronger 
(higher mass) than our final models described above and 
thus inconsistent with all of our multiple image systems. 
They also predicted an additional multiple image to the 
North which is not observed (in the vicinity of 17:22:27.8, 
+32:08:16). 

4.3. Complementary Strong Lensing Analyses 

We performed semi-independent lens modeling analy- 
ses using Lenstool (Kneib 1993; Julio et al. 2007) and 
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Figure 5. Projected 2D mass profiles from strong and weak 
lensing. We present seven SL models and two WL models. Five 
of the SL models are based on the Zitrin method (§4.1). One of 
these ( "Zitrin G" ) uses a different (Gaussian) kernel to smooth the 
light distribution for use as the halo mass model. The LensPerfect 
model is steepest in integrated M(< _R), which translates to the 
shallowest density profile k(R), while the Lenstool model is about 
average (§4.3). 

LensPerfect (Coe et al. 2008, 2010). In the course of 
these analyses, we verified the multiple image systems 
and estimated their redshifts independently. 

Our Lenstool model consisted of an NFW halo 
(Navarro et al. 1996) and truncated PIEMD (pseudo- 
isothermal elliptical mass distribution) halos (Kassiola & 
Kovner 1993) for the 69 brightest cluster members, which 
were again identified photometrically but independently 
from the analysis in §4.1. 

For the LensPerfect analysis, we assumed a prior that 
the mass is densest near the center of the BCG and 
roughly decreases outward radially. Otherwise, it in- 
cludes no assumptions about light tracing mass. Other 
priors include overall smoothness and rough azimuthal 
symmetry. For details, see Coe et al. (2008, 2010). The 
best solution found, according to these criteria, is shown 
in Fig. 4. It perfectly reproduces the observed positions 
of all multiple images (to the accuracy with which they 
are input) . We note the solution is only well constrained 
within the white polygon which bounds the multiple im- 
ages. 

4.4. Results from the Ensemble of Models 

Integrated projected mass profiles from our 7 SL mass 
models are presented in Fig. 5. We adopt the average 
and scatter of these models as our final SL constraints. 
Some level of agreement is guaranteed by the fact that all 
models used the same input multiple images and photo- 
metric redshift information. These image identifications 
were verified independently in each analysis. The red- 
shifts were allowed to vary somewhat and were optimized 
independently by each model. 

The models converge most tightly on the projected 
mass contained within ~ 20", roughly as expected given 
the Einstein radii of the systems. In Fig. 6, we plot the 
Einstein radius fl^ as a function of background source 
redshift z s for all 7 SL models. We calculated these 
Re(z s ) as those circular radii centered on the highest 
mass peak (coincident with the BCG) which enclose an 
average projected density equal to the critical strong 
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Figure 6. Einstein radius Re as a function of source redshift 
z s measured by our various lens models (with mass profiles plot- 
ted in Fig. 5). Our results are much lower than the rough esti- 
mate of Re ~ 40" (z s = 1.5) from ground-based imaging quoted 
by Umetsu et al. (2009). Instead, our models range between 
R E = 20" ± 2" (z s = 1.5) and 23" ± 2" (z a = 2). Some level 
of agreement is guaranteed by the fact that all models used the 
same input multiple images and photometric redshift information, 
though the rcdshifts are allowed to vary somewhat and are opti- 
mized independently by each method. The LensPerfect mass pro- 
file is shallow, hovering near the lensing critical density, allowing 
unconstrained perturbations beyond R > 30" (including some area 
outside the multiple image constraints) to significantly influence 
Re at these radii. 

lensing density. 

Our results, which range from Re = 20" ±2" (z s — 1.5) 
and 23" ± 2" (z s — 2), are significantly lower than those 
roughly estimated from ground-based imaging as quoted 
by Umetsu et al. (2009): R E = 40" ±4" (z s = 1.5). This 
previous estimate was based on Re ^ 30" for z s < 1 
assuming that the bright, prominent arc (marked with an 
X in Fig. 4) lies near the Einstein radius. In this work, 
we find this arc is not in fact located at R E - Instead, 
it lies at R m 27", greater than the R E = 17" ± 3" 
(z s = 1.2) determined robustly by our 12 other multiple 
image systems. Furthermore, we find it has a photo- 
z ~ 1.2, greater than that assumed in Umetsu et al. 
(2009). Both of these factors contributed to the higher 
concentration measured in that work. 

5. WEAK LENSING MASS MODELING 

To probe the mass distribution of A2261 at larger radii, 
we turn to weak lensing analyses of wider ground-based 
images obtained with Subaru Suprime-Cam, as previ- 
ously studied in detail by Umetsu et al. (2009) and Okabe 
et al. (2010). Here we present new, more robust analyses 
incorporating additional observations, improved image 
reductions, and new analytical techniques. The addi- 
tional observations include KPNO Mayall 4-m imaging 
(§5.1) and spectroscopy from MMT/Hectospec (§3). 

5.1. Subaru and KPNO Wide-Field Imaging 

Our weak lensing analysis is based on archival BjVjRc 
imaging obtained with the 34' x 27' FOV Suprime-Cam 
(Miyazaki et al. 2002) on the Subaru 8.2-m telescope, and 
NOAO archival MOSAIC1 i'z' imaging obtained with 
the Mayall 4-m telescope at KPNO (program 2008A- 
0356, P.I. Mandelbaum). The integration times are 20, 
30, 45, 60, and 90 minutes, respectively, for the five fil- 
ters BjVjRci' z' . The Subaru i?c-band images used for 



galaxy shape analyses have a seeing FWHM sa 0.65" and 
were obtained at two different orientations rotated by 90 
degrees (see §5.4). 

5.2. Previous Analyses 

Umetsu et al. (2009) and Okabe et al. (2010) both mea- 
sured weak lensing distortions in the Subaru i?c~band 
imaging described above. Background galaxies were se- 
lected based on VjRc color-magnitude cuts determined 
in part so as to minimize contamination of unlensed clus- 
ter galaxies as described in Medezinski et al. (2007) and 
Umetsu & Broadhurst (2008). 

Both analyses identified a massive background struc- 
ture at z ~ 0.5 which may affect the lensing signal. 
Umetsu et al. (2009) found that out of four clusters ana- 
lyzed, A2261 was the most sensitive to the exact profile 
fitting method used. 

The first method, used in both papers (and commonly 
elsewhere), fits the observed shears (binned radially) di- 
rectly to those expected from NFW profiles. The sec- 
ond method attempts to correct for the mass-sheet de- 
generacy based on the observed shears alone. With the 
outer "mass sheet" density left as a free parameter «(,, 
the discretized density profile is iteratively refined toward 
consistency with the observed reduced tangential shears 
g + = 7+/(l — k). This method (Umetsu & Broadhurst 
2008; Umetsu et al. 2009) is a non-linear extension of 
earlier "aperture densitometry" techniques developed by 
Fahlman et al. (1994) and Clowe et al. (2000). 

Umetsu et al. (2009) found that NFW fits to the pro- 
file derived from the latter method yield a marginally 
higher mass concentration c V i r = 10. 2j^'g than the 
c mr = 6A±\% obtained using the former method (Ta- 
ble 4 and Fig. 9). This higher value was supported by 
their rough estimate of the strong lensing Einstein radius 
Re ~ 40" (z s ~ 1.5) based on the prominent arc (see 
§4.4). A joint fit to the WL and R E = (40 ± 4)" yielded 
c vir = U.llf;!, which was quoted by Oguri et al. (2009) 
as an example of a cluster with higher than expected 
concentration. 

We revisit the Umetsu et al. (2009) c V i r ~ 10 WL result 
as follows. In that analysis, the innermost mass density 
profile point S(< 1') (corresponding to the mean den- 
sity interior to the inner radial boundary of weak lensing 
measurements) was not included in the fitting. Here we 
included that data point and found the best fit value 
and uncertainties both decreased from c V i r — IO.2I3 5 to 

c v i r = 5. 8^4. We note the significant improvements in 
precision gained by increasing the radial range of con- 
straints. 

Okabe et al. (2010) found a lower virial mass M V i r = 
1.36t^24 x 10 15 M Q and c vir = 6.04^;^. This profile un- 
derestimates mass in our strong lensing region by ~ 20% 
(§6). We compared our shear catalog directly with that 
used in Okabe et al. (2010) and provided to us. We find 
similar WL signal for R > 3' but recover stronger signal 
interior to this radius. This difference is most likely ex- 
plained by improved background selection in our catalog 
with lower contamination due to cluster galaxies (§5.5). 

5.3. Current Analysis 
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In addition to the Subaru Vj and i?c images, we also 
utilize the -Bj-band image, which improves our selection 
of background galaxies (§5.5) with respect to the pre- 
vious analyses. Our Subaru image reduction procedure 
(Nonino et al. 2009) is somewhat improved compared to 
that used in Umetsu et al. (2009) in terms of distortion 
corrections and image co-addition (here PSF- weighted). 
After trimming the shallower edges, the final co-added 
images roughly cover a circular area with a 17.8' radius 
(~ 1,000 square arcmin) which we use for our analysis. 
As in the previous analyses we measured galaxy shapes 
in the Subaru i?c-band images, though our procedure is 
slightly different (§5.4). 

The KPNO i' and z' images were reduced using calibra- 
tion frames, including fringe and pupil maps, obtained 
from the archive, and then stacked in a manner similar 
to the Subaru images. Zeropoints were calibrated based 
on comparisons with point source photometry from SDSS 
DR7 (Abazajian et al. 2009). 

Five-band BjVjRci'z' photometry was measured us- 
ing SExtractor in PSF-matched images created by Color- 
Pro (Coe et al. 2006). Subaru zeropoints were calibrated 
based on comparisons to HST and KPNO photometry, 
then recalibrated based on SED fits to photometry of 
galaxies with measured spectroscopic redshifts primarily 
from Hectospec and supplemented by SDSS DR7 (§3). 

5.4. Shape Measurement 

We produced two separate co-added i?c-band images 
for the shape analyses based on the imaging obtained 
at two different orientations separated by 90 degrees. 
Galaxy shapes (reduced shears) were measured at each 
orientation and their weighted averages computed. In 
both Umetsu et al. (2009) and Okabe et al. (2010), shapes 
were instead measured in co-added images which com- 
bined both orientations. We find this does not have a sig- 
nificant effect on the derived mass profile. When we use 
the same method to analyze shears measured in Umetsu 
et al. (2009) and in this work, we find consistent results 
(Table 4, rows 5 and 6). 

For accurate shape measurements of faint background 
galaxies, we used the IMCAT software (Kaiser et al. 
1995), following the formalism outlined in that paper. 
Full details of our weak lensing analysis pipeline are pro- 
vided in Umetsu et al. (2010). We have tested our shape 
measurement and object selection pipeline using simu- 
lated Subaru Suprime-Cam images (M. Oguri 2010, pri- 
vate communication; Massey et al. 2007) . We recover in- 
put WL signals with good precision: typically a shear cal- 
ibration bias \m\ < 5% (where this bias shows a modest 
dependence of calibration accuracy on seeing conditions) , 
and a residual shear offset c ~ 10 -3 , which is about one 
order of magnitude smaller than the typical distortion 
signal (reduced shear \g\ ~ 10 -2 ) in cluster outskirts. 
This level of performance is comparable to other similarly 
well-tested methods (Heymans et al. 2006), and has been 
improved in comparison with our previous pipeline used 
in Umetsu et al. (2009) which achieved 5% < \m\ < 10%. 

5.5. Background Galaxy Selection 

Robust selection of background galaxies is crucial in 
weak lensing analyses to minimize contamination by un- 
lensed cluster and/or foreground galaxies which would 
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Figure 7. Background galaxies are selected for weak lensing 
analysis (lower outlined region) based on Subaru BjVjRc color- 
color-magnitude selection. At small radius, the cluster overdensity 
is identified as the dashed region. The background samples are well 
isolated from this region and satisfy other criteria as discussed in 
§5.5. This background selection successfully excludes all 189 cluster 
members (black) within the Subaru FOV as identified based on our 
velocity caustic analysis of Hectospec spectroscopic redshifts (§3). 
Cyan points are Rq < 26 galaxies, where stars have been identified 
and excluded. 

dilute the lensing signal by a fraction equal to the level 
of contamination (Broadhurst et al. 2005b; Medezinski 
et al. 2007). If not accounted for properly, contamina- 
tion can be especially significant at small clustercentric 
radii where cluster galaxies are relatively dense. Previ- 
ous analyses have demonstrated that color-color selection 
using three Subaru broadband filters delivers robust dis- 
crimination between cluster, foreground and background 
galaxies (Medezinski et al. 2010, 2011; Umetsu et al. 
2010, 2011b,a). 

Here we began by detecting objects within 17.8' (~ 3.8 
Mpc) of the BCG (the area deeply imaged by Subaru). 
We pruned stars from this sample based on i?c-band 
magnitude, peak flux, FWHM, and SExtractor "stellar- 
ity". 

We then derived BjVjRq color-color- magnitude cuts 
(Figure 7) as described in Medezinski et al. (2007, 2010, 
2011). We calculated number count density and average 
clustercentric radius both as a function of position in 
this color-color space. Cluster galaxies are identified as 
a peak in the former and minimum in the latter. We de- 
termined the region occupied by these galaxies and later 
found it to coincide well with colors of cluster members 
as determined based on a velocity caustic analysis of Hec- 
tospec spectroscopic redshifts (§3). 

We then defined regions in this color-color space well 
separated from the cluster galaxies for use as our back- 
ground galaxy selection. The border placement is opti- 
mized to maximize total number counts while minimizing 
contamination from cluster members. The latter can be 
detected as dilution of the average shear signal and/or 
a rise in number counts toward the cluster center. We 
also imposed magnitude cuts 22 < Rq < 26 to further 
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Figure 8. Weak lensing tangential reduced shear (top) and mag- 
nification (number count depletion) profiles (bottom) measured in 
Subaru images of A2261 (§5.6). Also plotted is a joint Bayesian 
fit to both. This is our primary WL model. Bright objects are 
masked out to refine the estimates of area and thus number count 
densities. 

avoid contamination at the bright end and incomplete- 
ness at the faint end, while maximizing the number of 
faint galaxies which contribute to the lensing signal. 

Our final cuts (Fig. 7) yielded 12,762 background 
galaxies (12.8 / square arcmin) for WL analysis. We 
verified that the WL shear signal increases toward the 
center and that the B-mode (curl component) is consis- 
tent with zero. We later verified that these cuts success- 
fully reject all 189 galaxies identified spectroscopically as 
cluster members within the Subaru FOV (§3). 

To estimate the mean effective redshift of this back- 
ground population, we applied this same color-color- 
magnitude cut to galaxies with robust photometry and 
photometric redshifts measured in the COSMOS field 
(Capak et al. 2007; Ilbert et al. 2009). We compute the 
average lensing efficiency ft = D^s/Ds for this sample 
given our lens redshift zl = 0.225 and find an effective 
zs = 0.99±0.10 for our background WL sample. For each 
lensed galaxy, the factor ft is a function of angular diam- 
eter distances from lens to source D^s = Da( z Li z s) an d 
observer to source D$ = Da(0, zs)< We later marginal- 
ized over this uncertainty when fitting mass profiles to 
our WL data. 

5.6. Revised Weak Lensing Mass Profile 

In addition to our revised selection of background 
galaxies, we also used a slightly different method to esti- 
mate the "mass-sheet", or background density k^. Here 
we performed iterative NFW fitting allowing Kb to be a 
free parameter (Umetsu et al. 2010). 

We then performed a second analysis method which in- 
corporates independent weak-lensing magnification data 
(depleted number counts of faint background galaxies) 
in a Bayesian approach. For details on this method, 
see Umetsu et al. (201 lb, a). Measurements of number 
count depletion generally break the mass sheet degener- 
acy more robustly (e.g., Broadhurst et al. 1995; Umetsu 
et al. 2011b) than the aperture densitometry technique 
described above. 

We find a consistent mass profile solution based on a 
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Figure 9. Projected mass density profile derived from strong and 
weak lensing analyses fit to NFW profiles from published analyses 
(gray) and this work (colors). The SL results are the average and 
scatter of our seven mass models (§4.4). 

joint Bayesian fit to both the observed shears and magni- 
fication as shown in Fig. 8. The total signal-to-noise ra- 
tio (S/N) in our tangential distortion profile is S/N ss 17 
(defined as in equation 9 of Okabe et al. 2010), whereas 
S/N « 20 in the joint mass profile from combined tan- 
gential reduced shear and magnification measurements 
(equation 38 of Umetsu & Broadhurst 2008). Thus, in 
addition to breaking the mass sheet degeneracy, the mag- 
nification measurements also increased the overall signif- 
icance by ~ 20% (cf. Table 5 of Umetsu et al. 2011a; 
Rozo & Schmidt 2010). 

Our BjVjRc color-color selection does not allow us to 
effectively discriminate between "blue" and "red" back- 
ground samples with properties similar to those derived 
from BjRqz' color-color selection (Medezinski et al. 
2010, 2011). Galaxies in the "blue" samples identified in 
these works have steep number count slopes, roughly can- 
celing out any number count depletion. Stronger mag- 
nification signals are measured in "red" BjRqz' samples 
with relatively flatter number counts. 

To investigate the effect this may have on our analy- 
sis, we explored the BjRqz' colors of the subset of our 
galaxies detected in shallower z'-band KPNO imaging. 
We found that the majority of our background sample 
corresponds to a "red" selection in BjRqz' ', as desired. 
We repeated our magnification analysis on this red subset 
and found no significant changes in our results except for 
somewhat larger uncertainties due to the lower number 
of galaxies. 

6. MASS PROFILE FROM JOINT STRONG + WEAK 
LENSING ANALYSIS 

For a mass concentration determination that is both 
precise and accurate, the inner mass profile must be si- 
multaneously constrained by strong and weak lensing 
analyses (e.g., Meneghetti et al. 2010). In §4 we de- 
rived seven SL mass profiles from which we calculated 
the average with uncertainties. Then in §5 we presented 
various WL analyses. Our final WL analysis including 
both shear and magnification (number count depletion) 
information is our most robust. The magnification data 
break the mass sheet degeneracy and increase our overall 
WL signal-to-noise. 
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Table 4 

Published NFW Fits to the Mass Profile of A2261 
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Note. — We also roughly identify and estimate the effects of background/foreground structures along the line of sight (§6.1). We 
find that correcting for these may lower M v i T by ~ 7% and increase c v i r by ~ 5%. 
a The mass sheet density ftj, in the outer annulus was estimated based on an NFW fit to the shears. 
b In this work, we find a lower R E (z s = 1.5) = (20 ± 2)". 

c Now including the innermost bin k(< 1'), the mean convergence interior to the weak lensing measurements. 
d Iterative NFW fitting is performed to find the best fitting mass sheet density n^. 




Virial Mass M vlr (1O 15 M ) 

Figure 10. Constraints on the virial mass and concentration of 
A2261 from lensing analyses. Contours are l-cr and 2-<r (A% 2 = 
2.30 and 6.17, or fti 68% and 95% confidence assuming Gaussian 
ellipsoidal uncertainties). The dashed black line is l-cr for weak 
lensing if marginalizing over one of the variables (Ax 2 = 1)- Best 
fit values are also indicated. We marginalized over the weak lensing 
background redshift uncertainty z s = 0.99 ± 0.10. 

Our strong and weak lensing data agree well in their 
region of overlap (Fig. 9). We perform joint NFW fitting 
to the SL mass enclosed M{< R) measured at 12 points 
5" < R < 1' (18 - 215 kpc) and the WL mass density 
k(R) measured in 11 bins with centers 40" < R < 14.2' 
(144 - 3,059 kpc). 

This yields a virial mass M mr = 2.2±0.2 x 1O 15 M /i 7 " o 1 
and concentration c V i r = 6.2 ± 0.3 with a significantly 
greater precision than that obtained by WL alone (c V i r = 
5.7^J'7). Confidence contours are plotted in Fig. 10 and 
the constraints are tabulated in Table 4. Our new results 
strongly disfavor the previous c V i r ~ 10 results (Umetsu 
et al. 2009). 

Our use of the spherical NFW profile enables the most 



direct comparisons with analyses of simulated halos fit 
exclusively to this profile (e.g., Klypin et al. 2011; Prada 
et al. 2011; Bhattacharya et al. 2011). Other mass pro- 
files, including the Einasto (1965) profile, have been 
shown to yield slightly better fits to simulated halos 
(Navarro et al. 2004; Merritt et al. 2005, 2006; Navarro 
et al. 2010). The choice of profile does not significantly 
affect the derived concentrations (Duffy et al. 2008; Gao 
et al. 2008; Reed et al. 2011). 

6.1. Effect of Background Structures 

Significant background structure was identified by 
Umetsu et al. (2009) and Okabc et al. (2010) and es- 
timated to be at z ~ 0.5 based on Vj — Rc galaxy colors 
in that region. It was posited that this structure could 
bias the derived weak lensing signals. Here we present 
a rough estimate of the effects of background structures 
on our derived mass and concentration. 

We identified mass peaks in a weak lensing mass model 
obtained using a linear Kaiser & Squires (1993) mass re- 
construction method with Gaussian smoothing (Fig. 11). 
We then estimated redshifts for the twelve peaks with 
nearby bright galaxies based on spectroscopic and pho- 
tometric redshift information. Hectospec spectroscopic 
redshifts were available for ten of the peaks. For the 
remaining galaxies, we used BPZ photometric redshifts 
derived using their BjVjRci' z' magnitudes. We note 
these achieved a good accuracy of ~ 3%(1 + z) for the 
~ 300 galaxies with BjVjRci'z' photometry and confi- 
dent spectroscopic redshifts. 

We identified six mass peaks coincident with bright 
galaxies in the background or foreground. We then elim- 
inated those peaks from our mass model by setting the 
overdensity of those regions equal to zero (Fig. 11) and 
rederived the mass profile as determined by weak lens- 
ing. Based on fitting of NFW profiles to our strong and 
weak lensing, we found that removal of these background 
structures lowered the virial mass M V i r by ~ 7% and in- 
creased the concentration c V i r by ~ 5%. (Fitting to WL 
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Figure 11. Weak lensing mass reconstruction of A2261 based on analysis of Subaru images. Left: Mass contours are overlaid on 
this Subaru SjVj-Rc color image 27' ~ 5.8 Mpc on a side. Mass peaks are tentatively identified as belonging either to the cluster or 
background/foreground structures based on nearby galaxies with spectroscopic redshifts from Hectospec or 5-band Bayesian photometric 
redshifts (labeled with question marks "?"). Magenta labels correspond to background galaxies at 0.25 < z < 0.28 with recession velocities 
7,500 - 16,500 km/s greater than the cluster. Red labels correspond to 0.40 < z < 0.53. Right: Mass peaks are somewhat aggressively 
set to zero where contributions from background/foreground structures are suspected. The dark contour corresponds to zero projected 
overdensity. This is used to estimate the effects of large scale structure on our mass profile. We note this linear Kaiser &; Squires (1993) 
reconstruction was smoothed with a Gaussian kernel and not used in our analysis aside from this estimate. 



alone yielded slightly higher ~ 10% effects.) We con- 
clude that background structures likely affect the mass 
and concentration measurements from joint SL+WL fit- 
ting at the 10% level or less. We made some attempt to 
maximize this effect by setting mass overdensities equal 
to zero (some overdensity should remain in these regions 
due to the cluster). However our analysis was not ex- 
treme either in the number or sizes of areas eliminated. 



7. 



TRIAXIALITY FROM JOINT LENSING 
ANALYSIS 



X-RAY 



Lensing analysis may yield higher mass estimates than 
X-ray analysis for either or both of the following rea- 
sons: 1) halo elongation and/or additional massive struc- 
tures along the line of sight boosting the lensing signal 
(Meneghetti et al. 2010; Newman et al. 2011; Morandi 
&: Limousin 2011); 2) non-thermal gas pressure support 
(primarily turbulent flows and/or bulk motions) deviat- 
ing from assumptions of hydrostatic equilibrium (Nagai 
et al. 2007; Lau et al. 2009; Kawaharada et al. 2010). 

In cosmological simulations, dark matter halos are gen- 
erally found to be prolate with typical axis ratios of ~ 
2:1. This elongation is generally found to decrease as 
a function of radius (Barnes & Efstathiou 1987; Warren 
et al. 1992; Jing & Suto 2002; Schulz et al. 2005; Lemze 
et al. 2011). This trend may be dampened by baryons 
which are more dominant at smaller radii and act to 
make halos more spherical due to their collisional nature 
(e.g., Kazantzidis et al. 2004). Halo elongations along the 
line of sight can bias both lensing strengths and cluster 
concentration measurements significantly high, such that 
the measured concentrations of a lensing-selected sample 
may be biased high by ~ 50-100% (Hcnnawi et al. 2007; 
Oguri & Blandford 2009; Meneghetti et al. 2010, 2011). 



Non-thermal pressures may account for ~ 15% of the 
total support against gravitational collapse, thus biasing 
low by that amount the mass derived when assuming 
hydrostatic equilibrium (Nagai et al. 2007; Lau et al. 
2009). In relaxed clusters, non-thermal pressure support 
is expected to increase with radius up to ~ 30-40% at the 
virial radius due to inflowing gas (Lau et al. 2009; Shaw 
et al. 2010; Cavalierc et al. 2011b). A possible minimum 
in the non-thermal pressure support at ~ 0.1r V i r has also 
been predicted (Molnar et al. 2010). 

Previous joint lensing + X-ray analyses have allowed 
for these factors as global constants (e.g., Morandi et al. 
2010; Newman et al. 2011; Morandi et al. 2011b; Morandi 
<fc Limousin 2011). Radial dependence of non-thermal 
pressure support was modeled by Morandi et al. (2011a). 
Here we consider radial variation of this quantity as well 
as elongation. We only consider elongation along the 
line of sight as our 2D lens mass modeling already allows 
for elongation and more general asymmetries within the 
plane of the sky. 

7.1. Chandra X-ray Observations and Analysis 

A2261 was observed by Chandra ACIS-I in programs 
#550 and #5007 (P.I. Van Speybroeck) to depths of 
9.0 and 24.3 ksec, respectively (Morandi et al. 2007; 
Maughan et al. 2008; Gilmour et al. 2009; Mantz et al. 
2010a). 

We reprocessed and filtered the X-ray events in the lat- 
ter observation in a standard manner using CIAO v4.3 
and CALDB version 4.4.6. Based on — 24, 000 net pho- 
ton counts (0.7-7.0 keV), we extracted X-ray spectra 
within 15 annuli in the range 6" < R < 3.1' (20 kpc 
< R < 650 kpc) centered on the X-ray peak which is co- 
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incident with the center of the BCG. There were roughly 
equal net counts per annulus. A matched extraction 
of events from a reprojected, filtered, deep background 
events file was used for the background spectrum. 

XMM observed Abell 2261 on 9 separate occasions be- 
tween 2003 and 2004 for ~12-13 ksec each. Each ob- 
servation was heavily contaminated by proton flares and 
deemed unsuitable for analysis. It is likely that these 
lower priority observations were scheduled during peri- 
ods of elevated particle backgrounds. 

We fit the Chandra spectra simultaneously by creat- 
ing models of hot gas in hydrostatic equilibrium in a 
dark matter NFW gravitational potential well using the 
JACO (Joint Analysis of Cluster Observations) software 
(Mahdavi et al. 2007). JACO allows for nuisance param- 
eters such as an X-ray point source (none was detected) 
and contributions from a galactic soft background (found 
to be negligible in this case). 

In Fig. 12, we plot our NFW fit to the total mass (gas 
+ dark matter) profile assuming a spherical halo and hy- 
drostatic equilibrium (HSE). We fit out to r = 3.1' rj 667 
kpc hj Q , or just beyond r25oo = 590 kpc h^ , correspond- 
ing to Af 2500 = (0.29±0.05) x 10 15 M®h^, which we de- 
rive along with an NFW concentration C2500 = 2.3 ± 0.9. 
This mass is ~ 35% lower than the mass we derive at 
that radius based on our lensing analysis. For reference, 
if extrapolated to the virial radius, this profile would 
correspond to M vir = (0.82 ± 0.14) x 10 15 M Q h^ with 
c V i r — 9.1 ± 3.0. As we show, this is in good agree- 
ment with Zhang et al. (2010) who also fit out to larger 
radii using the XMM data. Maughan et al. (2008) find a 
slightly larger M 500 ~ 0.80 x 1O 15 M within R 50Q m 1.31 
Mpc based on the Chandra data. 

We also plot 20% deviations from HSE in the form of 
non-thermal pressure support. Though larger than ex- 
pected within r25oo (Lau et al. 2009; Shaw et al. 2010; 
Molnar et al. 2010; Cavaliere et al. 2011b; Nelson et al. 
2011), this is what the data would require to bring the 
lensing and X-ray masses derived in this work into agree- 
ment just within the error bars. 

Mantz et al. (2010a) derive a higher gas mass than 
Zhang et al. (2010) (also shown in Fig. 12). Based on 
this, they derive a significantly higher M500 = 1.44 ± 
0.26 x 10 15 M Q within r 500 = 1.59 ±0.09 Mpc. This is in 
excellent agreement with our derived lensing mass. 

Mantz et al. (2010a) assume a gas mass fraction f gas ~ 
12% for A2261, very similar to the f gas derived by Zhang 
et al. (2010) assuming hydrostatic equilibrium. Various 
systematics are discussed further in Conte et al. (2011) 
who also derive a range of mass estimates for A2261 sim- 
ilar to that described already. We consider this full range 
in our analysis. 

We note that published dynamical mass estimates of 
A2261 are significantly lower (Rines et al. 2010). These 
data are somewhat limited by bright stars in this field, 
hindering our ability to obtain additional spectra which 
might resolve this discrepancy. 

X-ray observables, and the masses derived from them, 
are largely insensitive to halo elongation (e.g., Gavazzi 
2005; Nagai et al. 2007; Buote & Humphrey 2011a,b). 
This is not the case for masses derived from lensing data, 
as we discuss below. 
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Figure 12. Mass profiles derived from various observational 
probes. We derive an X-ray mass (red curve, NFW profile) ~ 35% 
lower than our lensing mass (blue, NFW profile) at r2500 ~ 600 
kpc. Marginal agreement can be achieved by allowing for 20% non- 
thermal pressure support (magenta), though this is much higher 
than generally expected at this radius (e.g., Lau et al. 2009). Agree- 
ment may be more readily achieved by an elongated halo with a 2:1 
axis ratio for r > 100 kpc (light blue hashed). However, the need 
for such elongation may be obviated completely by systematic un- 
certainties in the X-ray results (Maughan et al. 2008; Zhang et al. 
2010; Mantz et al. 2010a). The latter result agrees well with our 
spherical lensing mass profile. A similar range of results was found 
by Conte et al. (2011) who consider various systematics. The dy- 
namical mass estimates Rines et al. (2010) are significantly lower 

(Mioo ~ 0.5 - 0.7 x lO 1B M Af o 1 ). Also plotted are gas mass mea- 
surements based on X-ray (orange) and AMiBA SZE observations 
(brown stars; Umetsu et al. 2009). 

7.2. Halo Elongation 

As shown in Fig. 12, our mass profiles derived inde- 
pendently from lensing and X-ray analyses are in good 
agreement in the core, while the latter exhibits a ~ 35% 
deficit at the X-ray r 2 soo ~ 600 kpc. This result is to- 
ward the low end of other X-ray mass estimates, so we 
consider this to be a limiting case. This deficit could best 
be accounted for by halo elongation along our line of sight 
(though the Mantz et al. 2010a result would obviate the 
need for any such elongation). 

We found that an axis ratio of 2:1 is able to bring our 
lensing and X-ray results into better agreement at T25oo- 
This elongation is not required at inner radii where a 
spherical profile fits the data. Halo elongation is gener- 
ally expected to decrease, not increase, with radius (e.g., 
Hayashi et al. 2007). However we note that here we 
are probing the very inner core where the dense con- 
centration of baryons may increase the sphericity (e.g., 
Kazantzidis et al. 2004). The large BCG of A2261 ex- 
tends visibly to r ~ 100 kpc. 

We construct a toy model for the halo elongation e = 
1 — b/a varying with radius, increasing from zero (spher- 
ical) for r < 1 kpc to 0.5 (an axis ratio of 2:1) beyond 
r > 100 kpc. Between these two radii, it follows e(r) = 
0.251og 10 (r/kpc). The 3D mass density p(r) scales with 
halo roundness (£ = 1 — e): p{r) = £(r)/?NF\v(u) where 
u = \lx 2 + y 2 + £z 2 and pnfw(u) = /9 s (u/r s ) _1 (l + 
u/r s )~ 2 . This scaling preserves the projected mass den- 
sity k(R) integrated along the line of sight (z-axis) and 
thus preserves all lensing observables. 
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We applied this elongation profile to our primary joint 
(SL + WL shear + magnification) lensing profile NFW 
fit: M vir = 2.2 ± 0.2 x W 15 M Q h^ and c vir = 6.2 ± 0.3. 
We then calculated numerically the 3D mass enclosed 
within spherical shells for this ellipsoidal mass distribu- 
tion. This is plotted as the light blue hashed region in 
Fig. 12. We find this model agrees well with both our 
lensing and X-ray derived mass profiles, whether includ- 
ing modest non-thermal pressure support or not. 

A spherical NFW fit to this elongated profile yields 
M mr = 1.7 ±0.2 x lO 15 M hm and c vlr = 4.6 ±0.2. Ap- 
plying the corrections for background / foreground line 
of sight structures estimated in §6.1, we find M V i r ~ 
1.6 x 10 15 Af Q (a - 7% decrease) and c vir ~ 4.8 (a - 5% 
increase). Note the former corrections for cluster halo 
elongation are significantly larger than those for LOS 
structure. 

We note that this measurement method is consistent 
with that generally used to measure the mass profiles of 
simulated clusters. Enclosed mass (or, more often, den- 
sity) is determined assuming spherical symmetry (and 
most often fit to an NFW profile) even though the halos 
are triaxial and asymmetric. Another approach is to fit 
the lensing and X-ray observables to an ellipsoidal NFW 
profile, as in the Morandi & Limousin (2011) analysis of 
Abell 383, the first observed CLASH cluster. Notably, 
they allow for a fully general ellipsoidal gNFW (gener- 
alized NFW with variable inner slope) dark matter halo 
plus an exponential ICM profile including non-thermal 
pressure support. Ideally, simulations will be analyzed 
in the same way allowing for direct comparisons. Un- 
til then, the advantages of this parameterization will 
not be completely realized, as spherical averages must 
be derived for comparison with most published analy- 
ses of simulations. Morandi & Limousin (2011) derive 
M mr = 8.6 ± 0.7 x 10 14 Af© and c vir = 6.0 ± 0.6 (pri- 
vate communication) based on a joint SL + X-ray anal- 
ysis. We compare this to the M vir = 7.7 ± 1.0(stat.) ± 
O.4(syst.)xlO 14 M and c mr = 8.8±0.4(stat.)±0.2(syst.) 
found by Zitrin ct al. (201 Id) who fit a spherical NFW 
profile to joint SL + WL constraints. The effect of cor- 
recting for elongation is to decrease the derived concen- 
tration as in our analysis of A2261 (see Fig. 13). 

8. MASS PROFILE COMPARED TO SIMULATED 
CLUSTERS 

Based on our joint strong + weak lensing + X-ray 
analysis, we find that A2261 is not significantly over- 
concentrated compared to simulated relaxed clusters of 
similar mass and redshift. This is demonstrated in 
Fig. 13. Our range of allowed concentrations (4.4 < 
c I^j 6.5) span the low end of average expectations 
(4-5 < (c) < 7.8) from simulations (Duffy ct al. 2008; 
Zhao et al. 2009; Klypin et al. 2011; Prada et al. 2011; 
Bhattacharya et al. 2011). Results are also plotted for 
the first CLASH cluster Abell 383 (Zitrin et al. 2011d; 
Morandi & Limousin 2011). 

Note the recent Bolshoi and Multidark simulations 
(Klypin et al. 2011; Prada et al. 2011) produce halos 
with significantly higher concentrations than previous 
simulations (Neto et al. 2007; Maccio et al. 2008; Duffy 
ct al. 2008; Zhao et al. 2009, although see Hcnnawi et al. 
2007). Prada ct al. (2011) find upturns in c{M,z) for 
high masses and redshifts. This behavior is unexpected 
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Figure 13. Observed virial masses M v i T and concentrations c v i r 
for CLASH clusters compared to the average c(M, z) realized for 
relaxed clusters in simulations. Squares are from joint strong + 
weak lensing analyses of A2261 (this work) and A383 (Zitrin et al. 
2011d). The hexagon is from Morandi & Limousin (2011, and 
private communication) who fit triaxial halos to A383 SL + X- 
ray data. For A2261, we plot both error bars (1-cr, marginalizing 
over the other parameter) and confidence contours (l-<r). System- 
atic uncertainties are labeled: possible halo elongation (§7.2) and 
line of sight structures (§6.1). Results realized in two simulations 
(Duffy et al. 2008; Prada et al. 2011) are shown in light blue, in- 
cluding scatters of ~ 0.1 in log 10 (c) (~ 26%). Portions of these 
lines are dashed to indicate extrapolations to high masses where 
clusters are not realized in sufficient numbers. Averages results 
from three additional simulations (Zhao et al. 2009; Klypin et al. 
2011; Bhattacharya et al. 2011) are shown in red with styles solely 
for clarity. Results are plotted for relaxed cluster subsamples as 
determined by Duffy et al. (2008) and Bhattacharya et al. (2011), 
yielding concentrations ~ 10% higher than for the full populations. 
This 10% factor is applied to the results from the other simulations. 
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Figure 14. Similar to Fig. 13 but now for all clusters (not just re- 
laxed). The additional data points are all analyses of non-CLASH, 
lensing-biased clusters, as follows and all colored accoring to red- 
shift. Circles are from Umetsu et al. (2011b) SL+WL analyses. 
And small diamonds are from Oguri et al. (2011) analyses with 
WL + basic SL constraints (the Einstein radii). Note the average 
predictions from simulations for c(M, z) for all clusters are ~ 10% 
lower than for relaxed clusters. The expected scatters are larger: 
~ 0.15 in log 10 (c) (~ 41%). 
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and its origin needs to be understood. Bhattacharya 
et al. (2011) find no evidence for such upturns in their 
analysis of an even larger simulation, 2 Gpc h^ 1 on a 
side, 8 times the volume of Multidark, with the same 
number of particles (2048 3 ). 

A383 and A2261 are two of 20 CLASH clusters selected 
based on X-ray properties. We expect this sample to be 
less biased toward elongations along the line of sight than 
a lensing-selected sample. However some bias may re- 
main. If we assume all clusters are prolate and elongated 
in some direction, perhaps clusters which are roughly 
round and symmetric in the X-rays may preferentially 
be elongated along our line of sight. We will continue 
to make precision measurements of the mass profiles and 
constrain DM elongation for the CLASH clusters as the 
survey progresses. 

A2261 is borderline relaxed. Gilmour et al. (2009) clas- 
sified it as disturbed, but the X-ray peak is well aligned 
with the BCG, and the X-ray luminosity is symmet- 
ric except for a subclump to the SW. Maughan et al. 
(2008) measured centroids in various annuli and found 
the RMS shift to be modest (w) = (7.1 ±0.6) x 10- 3 i? 50 o, 
consistent with that found for relaxed clusters (w) < 
0.012i? 500 . 

In Fig. 14, we plot the expected c(M, z) for all clusters 
versus the most robust results from other strong + weak 
lensing analyses to date, including those just mentioned 
plus Umetsu et al. (2011b) and Oguri et al. (2011). These 
clusters were initially selected for study based on their 
lensing strength, thus their concentrations are expected 
to be biased significantly high. Disparity in this compar- 
ison is further increased as the expectations from sim- 
ulations are lower. Average concentrations for relaxed 
clusters (as plotted in the previous figure 13) are found 
to be ~ 10% higher (and have lower scatter) than aver- 
ages for the general population as plotted in this figure 
14. 

9. DISCUSSION: THE FORMATION TIME OF A2261 

We have found the mass profile and concentration of 
A2261 to be in agreement with values realized in cosmo- 
logical simulations for similar clusters. This is contrary 
to the previous finding of Umetsu et al. (2009) based 
solely on ground-based data which found a high concen- 
tration suggesting an early formation time. 

We can attempt to quantify this statement based on 
the relation c b s w c\a b s jaf found in previous work 
(Bullock et al. 2001; Wechsler et al. 2002; Zhao et al. 
2003; Wechsler et al. 2006). Here a = (1 + z)- 1 is the 
cosmic scale factor. The halo concentration is imprinted 
with the background density at its formation time via 
c(zf) w ci(l + Zf) and then increases over time roughly 
as c{z) oc a = (1 + z) -1 . 

The constant ci depends on the criteria used to de- 
fine "formation time" . Various definitions have been pro- 
posed based on the slowing rate of mass accretion (e.g., 
Wechsler et al. 2002; Cavaliere et al. 2011a) or mass at- 
taining some fraction of the observed mass (e.g., Sadeh 
& Rephaeli 2008). 

We present results based on these various definitions 
in Table 5. Regardless of the definition, we note that 
the Umetsu et al. (2009) result of c V i r ~ 11 implies a 
formation time (1.7 < Zf < 2.9) several billion years 



Table 5 

A2261 Formation Redshift zj Estimates 
Based on Various Criteria for "Formation" 
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a Umetsu et al. (2009) result. 
b This work: spherical halo. 
c This work: elongated halo. 

d Normalization according to each definition, respec- 
tively Cavaliere et al. (2011a); Wechsler et al. (2002); 
Sadeh & Rephaeli (2008). 

earlier than our primary result for a spherical halo c V i r ~ 
6.3 (0.5 < Zf < 1.2). This, in turn, implies a formation 
time several billion years earlier than our result for an 
elongated halo c V i r ~ 4.6 (0.13 < Zf < 0.6). The lone 
Zf < z b s — 0.225 result would suggest the cluster has yet 
to finish "forming" according to the Sadeh & Rephaeli 
(2008) definition. 

Concentration may be the observable most tightly cor- 
related with age for relaxed clusters (Wong & Taylor 
2011), but other probes may also be brought to bear. 
Smith et al. (2010) studied BCG morphology, luminosity 
gap Ami2 between the brightest and second brightest 
cluster galaxy, substructure fraction f su b, and cool core 
strength, as well as concentrations (as available from X- 
ray analyses in Sanderson et al. 2009) in a sample of 59 
massive clusters, including A2261. A2261 was found to 
be one of four "fossil clusters" with a large luminosity gap 
Ami2 > 2. Clusters with Ami2 > 1 were found to have 
less substructure, stronger cool cores, and higher mass 
concentrations, all likely signatures of earlier formation 
times without recent major mergers. In these clusters, 
the BCG has presumably had time to grow and accrete 
a significant fraction of the substructure mass (see also 
Ascaso et al. 2011). 

Based on X-ray observations, A2261 is borderline re- 
laxed (see discussion in §8) and a borderline cool core 
cluster. Though the temperature profile dips down in 
the core (Cavagnolo et al. 2009), the central entropy floor 
(Kq — 61 ± 8 keV cm 2 ) is higher and the density profile 
slope (a ~ —0.7 at 0.04 r^oo) shallower than generally 
found (Kq < 30 keV cm 2 and a < —0.85) for cool core 
clusters (Sanderson et al. 2009; Smith et al. 2010). There 
is no obvious star formation visible in the NUV/optical as 
often found in cool core clusters. A radio source aligned 
with the BCG is detected with - 5.3 mJy ~ 8 x 10 23 
W Hz" 1 at 1.4 GHz in NVSS (Condon et al. 1998) and 
3.39 mJy at 21 cm in FIRST (Becker et al. 1995). All 
69 radio-bright (> 2 x 10 23 W Hz" 1 at 1.4 GHz) BCGs 
analyzed by Sun (2009) were found to be in X-ray cool 
cores. 

Ultimately, analyses of these various observables in 
CLASH clusters and in simulated clusters with similar 
properties will contribute to significant advancements in 
our understanding of structure formation and evolution. 

10. CONCLUSIONS 



CLASH: Precise Mass Profile of Abell 2261 



We performed the first robust joint strong and weak 
lensing analysis of the galaxy cluster A2261. We find 
a halo virial mass M vir = 2.2 ± 0.2 x 10 15 M & h^Q and 
concentration c V i r = 6.2 ± 0.3 when assuming a spheri- 
cal halo. These tight constraints were enabled through 
a combination of the 16-band imaging from CLASH 
with multiband wide-field imaging from the Subaru and 
KPNO telescopes. The results show that A2261 is not 
"over-concentrated" as previously found but rather is in 
good agreement with predictions from ACDM N-body 
simulations. 

To explore halo elongation along the line of sight, we 
also derived a mass profile based on Chandra X-ray data, 
finding it to be ~ 35% below the lensing mass profile at 
^2500 (~ 500 kpc). This deficit may be explained by an 
axis ratio of ~ 2:1 outside the inner core r ~ 100 kpc, cor- 
responding to the visible extent of the BCG. This elon- 
gated mass profile has a lower spherically-defined virial 
mass M V i r = 1.65i{j!i2 x W 15 M & h^ and concentration 
c V ir — 4.6±0.2. Correcting for the lensing effects of mas- 
sive background structures may increase c V i r by ~ 5% 
and decrease M vir by ~ 7%. This lower c„j r ~ 4.8 
still agrees with predictions from many simulations but 
is lower than predicted by one recent study (Prada et al. 
2011). 

The need to assume halo elongation is critically tied 
to the reliability of the X-ray mass profile. Non-thermal 
pressure support may account partially for the lower X- 
ray mass. Published X-ray mass estimates have signifi- 
cant scatter, including one result in excellent agreement 
with our spherical lensing mass at r§QQ (^1.6 Mpc). 

The CLASH survey is providing fundamental and sub- 
stantial improvements in the quantity and quality of 
observational constraints on cluster dark matter halos. 
Simulations will be tasked with reproducing these em- 
pirical results, contributing significantly to our under- 
standing of structure formation. Ultimately our results 
will either confirm ACDM predictions or perhaps yield 
clues as to the nature of dark energy. 
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